Geometric phases and quantum phase transitions in open systems 
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tum system governed by the non-Hermitian effective Hamiltonian with the accidental crossing of 
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demonstrated that related quantum phase transition is of the first order. 
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Quantum phase transition (QPT) is characterized by 
qualitative changes of the ground state of many body 
system and occur at the zero temperature. QPT being 
purely quantum phenomena driven by quantum fluctu- 
ations is associated with the energy level crossing and 
implies the lost of analyticity in the energy spectrum at 
the critical points [1]. A first order QPT is determined 
by a discontinuity in the first derivative of the ground 
state energy. A second order QPT means that the first 
derivative is continuous, while the second derivative has 
either a finite discontinuity or divergence at the criti- 
cal point. Since QPT is accomplished by changing some 
parameter in the Hamiltonian of the system, but not the 
temperature, its description in the standard framework of 
the Landau- Ginzburg theory of phase transitions failed, 
and identification of an order parameter is still an open 
problem [2]. In this connection, an issue of a great in- 
terest is recently established relationship betweengeo- 
metric phases and quantum phase transitions @> 0> HQ- 
This relation is expected since the geometric phase as- 
sociated with the energy levels crossings has a peculiar 
behavior near the degeneracy point. It is supposed that 
the geometric phase, being a measure of the curvature of 
the Hilbert space, is able to capture drastic changes in 
the properties of the ground states in presence of QPT 

mm. 

In this Rapid Communication we analyze relation be- 
tween the geometric phase and QPT in an open quantum 
system governed by non-Hermitian Hamiltonian. We 
found that QPT is closely connected with the geometric 
phase and the latter may be considered as an universal 
order parameter for description of QPT. Studying the 
dissipative one-dimensional Ising model in a transverse 
magnetic field we demonstrated that the QPT being of 
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the second order in absence of dissipation is of the first 
order QPT for the open system. 

Degeneracy points and geometric phase.- We consider 
an open quantum mechanical system which together with 
its environment forms a closed system. The description 
of the such systems by effective non-Hermitian Hamilto- 
nian is well known beginning with the classical paper s by 
Weisskopf and Wigner on the metastable states [8|, [9] [30] . 

For the Hermitian Hamiltonian coalescence of eigenval- 
ues results in different eigenvectors, and related degen- 
eracy referred to as 'conical intersection' is known also 
as 'diabolic point'. However, in a quantum mechanical 
system governed by non-Hermitian Hamiltonian not only 
merging of eigenvalues of the Hamiltonian but the asso- 
ciated eigenvectors can be occurred as well. The point of 
coalescing is called an "exceptional point". At the lat- 
ter the eigenvectors merge formi ng a Jordan block (for 
review and references see e.g. [Tol flTj). 

In the context of the Berry phase the diabolic point 
is associated with 'fictitious magnetic monopole' as fol- 
lows. Assume that for adiabatic driving quantum system 
two energy levels may cross. Then the energy surfaces 
form the sheets of a double cone, and its apex is called a 
"diabolic point" [12[ . Since for generic Hermitian Hamil- 
tonian the codimension of the diabolic point is three, it 
can be characterized by three parameters R = (X, Y, Z). 
The eigenstates |n, R) give rise to the Berry's connec- 
tion defined by A n (R) = z(n, R|Vr|ti, R), and the cur- 
vature B n = Vr x A n associated with A n is the field 
strength of 'magnetic' monopole located at the diabolic 
point pjj 03 • The Berry phase j n = § c A n • dH is 
interpreted as a holonomy associated with the parallel 
transport along a circuit C [15]. Similar treatment of the 
non-Hermitian Hamiltonian yields the 'fictitious complex 
monopole' located at the exceptional point [16| . 

For the first time, the extension of the Berry phase 
to the non-Hermitian systems has been done by Gar- 
rison and Wright as follows [17]. Let an adjoint pair 
be a solution of the time dependent 
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Schrodinger equation and its adjoint equation (H = 1) 



-i^(*(t)\ = m)\H(\(t)), 



(i) 

(2) 



where A G 9Jt, the parameter space being 9JX. Let 
|^ n (A)) and (^ n (A)| being right and left eigenvectors 
of the Hamiltonian: iJ(A)|^ n (A)) = E n (X)\tp n (X)) , 
(ijj n (X)\H(X) — E n (X)(ip n (X)\. Now suppose that there 
exists a time period T for which A(T) = A(0), then a com- 
plex geometric phase j n is given by the integral [13, G3] 



7n 



A (n) 



(MX)\v a \M^)dx a 

(^n(A)|^n(A)) 



(3) 



where the integration is performed over the contour C in 
the parameter space, a = 1, . . . , dim9JT, being the 
connection one-form . Further we assume that the in- 
stantaneous eigenvectors form the bi-orthonormal basis, 

$m\^n)=Smn MM- 

Geometric phase and quantum phase transition. 
Analysis of the relation between QPT and geometric 
phase we begin with consideration of a two-level sys- 
tem described by generic non-Hermitian Hamiltonian 
H = Aol + R(£) • 0", where G{ are the Pauli matrices, 
R(t) = (X, y, Z) is slowly varying and Ao, X, F, Z G C. 
Using the spinless fermionic creation and annihilation 
operators, which obeys an anticommutation relations 
{C,C} = 0, {C+,C+} = 0, and {C,C+} = 1, one can 
rewrite the Hamiltonian as H = (Ao — i?)l + 2RC^C, 

where R = (X 2 + Y 2 + Z 2 ) 1 ' 2 . The ground state \u-) is 
defined as the vacuum state determined by C\u-) = 0. 
The instantaneous eigenvectors are found to be 



sin : 



COS : 



sin : 



. <5-| = 



e * sin - , cos - 
2' 2 



2 ,<s+l = 



cos-, sin- ] (4) 



where are the complex angles of the complex spher- 
ical coordinates, and the complex energy spectrum is 
given by E± = Ao =L i?. Coupling of eigenvalues occurs 
when R = and there are two cases. The first one is of 
the diabolic point located at the origin coordinates. The 
second case yields the exceptional point (Xq, Yq, Zq). At 
the latter the eigenvectors coincide up to the phase fac- 
tor, \u+) = e iK \u_) and (u + \ = e~ iK (u-\ 

The geometric phase of the ground state is given by 
7 = (1/2) § c q(l — cos6)d(p, where integration is per- 
formed over the contour C on the complex sphere S 2 . 
Let us assume that the contour C of integration is chosen 
as 6 = const. Then the geometric phase of the ground 
state is given by 7 = 7r(l — Z/R) and can be written 
as 7 = 7r(l + dE-/dZ), where E- is the ground state 
energy. As can be observed, lost of analyticity occurs at 
the degeneracy 'point' defined by R = and on the Dirac 



string attached to the complex fictitious monopole and 
crossing the complex sphere S 2 at the south pole. 

Further simplification can be made writing R = p — is, 
where we set p = (x, y, z). Without loss of generality we 
may choose the coordinate system thus, that e = (0, 0, e). 
Then computation of geometric phase yields 



7 = 7T 1 



z — is 



\]r 2 + (z — ie)* 2 



(5) 



where r = ^/ x 2 + y 2 . 

In what follows we consider the behavior of the geo- 
metric phase near the critical points, starting with the 
Hermitian Hamiltonian. Inserting e = in Eq. ([5]), we 
obtain 7 = 7r(l — z/(r 2 + z 2 ) 1 ! 2 ). This implies that the 
geometric phase behaves as the step-function near the 
diabolic point. Considering the general case, we obtain 



Re7 = 



7r, if r > £, (z = 0) 



K lT vP=3r)' ifr<£ '* 



±0 



(6) 



I e A — r* 

where the upper/lower sign corresponds to z — > ±0, As 





FIG. 1: Left panel (e = 0): There is clear step function be- 
havior of the geometric phase at the diabolic point r = z — 0. 
Right panel: Re7 nearby the exceptional point (e = 0.5). 

can be observed in Fig. [TJ if e = the geometric phase 
behaves as step-like function near the diabolic point. In 
addition, Re7 — > ±00 at the exceptional point r = e, and 
it behaves as a step-like function while r — > 0. Similar 
consideration of the imaginary part yields 



Ini7 



0, if r < e (z = 0) 

, o 2 > if r > e (z = 0) 

— £ 



(7) 



r 2 — e 2 

and clearly it diverges at the exceptional point, Inry 



Once return to general non-Hermitian X-dimensional 
problem, consider the non-Hermitian diagonalizable 
Hamiltonian H(X) = Yli=i ^I^X^I- The ground state 
is given by \i/j g (X)) = ^fL^ipi (A)), and computation of 
geometric phase yields 



7 = z^(^(A)| — |V a (A))dA a 



N 

£7* 



(8) 
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where 7^ is the geometric phase associated with the eigen- 
vector |^(A)). Then, applying the Stokes theorem and 
the Sghrodinger equation H\xm) = Em\Xm) together 
with its adjoint equation (x m \H = E m (xm\i we obtain 



N N 



(K\VaH\Xm)(x m \VbH\ Xl )d\ a A d\ b 



It follows herefrom that the curvature = <L4.W di- 
verges at the degeneracy points, where the energy levels, 
say E n and E n +i, are crossing, and reaches its maximum 
values at the avoided level crossing points. Thus, the crit- 
ical behavior of the system is reflected in the geometry 
of the Hilbert space through the geometric phase of the 
ground state. 

Since in the neighborhood of either diabolic or excep- 
tional point only terms related to the invariant subspace 
formed by the two-dimensional Jordan block make sub- 
stantial contributions, the TV-dimensional problem be- 
comes effectively two-dimensional (for details see [H, 
|20| ) . This implies that there exists the map (p : 9DT S 2 
such that in the vicinity of the degeneracy points the 
quantum system can be described by the effective two- 
dimensional Hamiltonian H e f = Aol + R • <x, where 
R = (E n +i — E n )/2. Then we have 
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R-dS 



7i(R) 



(9) 



where E' = y?(E) C S 2 . The behavior of the geomet- 
ric phase described by the first term is independent of a 
peculiarities of quantum-mechanical system. Therefore, 
one can consider the complex Bloch sphere as an univer- 
sal parameter space for description of QPT in the vicinity 
of the criticalpoint. 

Following Q, we define the overall geometric phase of 
the ground state as 7^ = ^ n ^ ne ther- 

modynamical limit j g = J j(x)dfi(x), where dfitx) is the 
suitable measure. As has been shown by Zhu [5[ on ex- 
ample of XY spin chain, the overall geometric phase as- 
sociated with the ground state exhibits universality, or 
scaling behavior in the vicinity of the critical point. In 
addition, the geometric phase allows to detect the crit- 
ical point in the parameter space of the Hamiltonian 
@, El HI S 0] • These works indicate that the overall geo- 
metric phase j g can be considered as the universal order 
parameter for description of QPT. 

Geometric phase and QPT in the quantum Ising model. 
- As illustrative example we consider the 1-dimensional 
Ising model in a transverse magnetic field with dissipa- 
tion governed by the non-Hermitian Hamiltonian: 

H = - J H ( k < + «+l ~ ^>n) (10) 
n=l Z 

with the periodic boundary condition cr^v+i = (T\. The 
external magnetic field is described by the parameter h 



and spontaneous decay is described by Y = V^cr~ with a 
source of decoherence being a~ = (cr* — zcr^)/2. 

To study the geometric phase in this system we con- 
sider the more general Hamiltonian H(h, J, ip) — g^Hg^ 

where = Yl^=i e % ^ a ™ and < cp < 2tt. After applying 
the standard Jordan- Wigner transformation and follow- 
ing the procedure outlined in 0, [2l|, we find that the 
system can be described in terms of the non-interacting 
quasiparticles with the reduced Hamiltonian 



N 



H + = - J ^2 ( c n c n+i + e l(p c n+1 c n +g + i5 



2gc\ l c n + c^ nJrl c n 



C n C n+l) 



(11) 



where g = h — i5, and c n are fermionic operators sat- 
isfying anticommutation relations {c m ,c^} = S mn and 
{c m ,c n } = {c^,^} = 0. Applying the Fourier trans- 
formations c n = e- in ^^2k c k eikna / Nl/2 with the an- 
tiperiodic boundary condition cat+i = — c±, we obtain 

H + = JY.k (% - cos(fca))4c fc + sin(fca)(e-^4c! fe + 



°C-kCk) 



z^l, where k = ±ir/Na, . . . ,±(N 



l)ir/Na is a half-integer quasimomentum, the lattice 
spacing being a. 

The Hamiltonian H + can be diagonalized by us- 
ing the Bogoliubov transformation: c& = UkH + 
v-k^_ k , c\ = u k b\ + v-kb-k- The Bogoliubov modes 
(iik,Vk) and (uk,Vk) satisfy the Schrodinger equation 
and its adjoint equation, respectively, with the Hamil- 
tonian H(k) = -iJSl + R(fc) • tr, and R(fc) = 
2 j( sm(ka) cos ip, sm(ka) sin cp, g — cos(fca)). There are 
two eigenstates for each k with the complex energies 
e±(k) =£o± where we set Sq = — iJ5 and e(k) = 
2J(g 2 — 2gcos(ka) + l) 1 / 2 . The positive energy eigen- 

state \u+(k)) = ^ ^ ^ > = normalized 

so that UkUk + VkVk = 1, defines the quasi-particle oper- 
ators b k = u k c k + VkcL k and b\ = u k c\ + v k C- k as fol- 
lows: b k = cos \ c k + sin \ cl k , b\ = e-** cos ^ c\ + 
sin ^ c_fc, where 



cos 9 k 



cos(ka) 



\/ g 2 — 2gcos(ka) + 1 



(12) 



Using these results, we obtain the diagonalized Hamil- 
tonian as a sum of quasi-particles with half-integer quasi- 
momenta, H + = ^ k (e + e(k)(b\b k - |)). Its ground 
state is given as product of qubit-like states: 

l^)=(g)(cos^|0} fc |0)_ fc -e-^sin^|l} fc |l)_ fc ) 

k 

$ g \ = (g) (cos^(0| fe (0|_ fe - e^sin^(l| fe (l|_ fe ) 

k 

where \0) k is the vacuum state of the mode bk, and \l)k 
is the first excited state, \l)k = b\,\0)k- Each single 
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state lies in the two-dimensional Hilbert space spanned 
by \0)k\0)-k an d For given value of k the state 

in each of these two-dimensional Hilbert space can be 
presented as the point on the complex two-dimensional 
sphere S% with coordinates (9k, <p)- 

For \g\ ^> 1 the ground state is a paramagnet with all 
spin oriented along the x axis, and from Eq. ([12]) we ob- 
tain cosOk — > 1 while \g\ — > oo. Thus, the north pole of 
the complex Bloch sphere corresponds to paramagnetic 
ground state. On the other hand, when \g\ <C 1 there 
are two degenerate ferromagnetic ground states with the 
all spins polarized up or down along the z axis. The real 
part of the complex energy reaches its minimum at the 
point defined by cos 9k = — 1, and, hence, the south pole 
of the complex sphere is related to the pure ferromagnetic 
ground state with the broken symmetry when all spins 
have orientation up or down. However, in the thermo- 
dynamical limit the system passing through the critical 
point ends in a superposition of the up and down states 
with finite domains of spins separated by kinks [2lj . 

The geometric phase of the ground state is found to be 

7 = */ $g\^W 9 ) dtp = ^n(l-cos0 k ) (13) 
Jo ^ k>0 

As can be shown, in the thermodynamical limit the 




FIG. 2: Real part of the overall geometric phase j g (left) 
and imaginary part of the overall geometric phase j g (right) 
versus S and ft. 




FIG. 3: Real part of derivative of the overall geometric phase 
djg/dh (left) and its imaginary part (right) versus 6 and ft. 

energy gap As (ft, k) vanishes and the geometric phase 



diverges at the exceptional point ft c = (1 — £ 2 ) 1//2 , 
k c = arcsin^/a. However, the overall geometric phase 
lg ~ ( n /N) Sfc>o(l — cos9k) written in thermodynami- 
cal limit as 

C I Q — cos x \ 7 , v 

^ = / 1 " 7=2 9 ^ (M) 

y v g z — Zgcosx + I J 

has finite jump discontinuity at the exceptional point 
(Fig. [2|). The result of integration can be written in 
terms of the complete elliptic integrals of the first and 
second kinds 




We note that j g can be written as j g = 7r(l + dE g /dh), 
where E g = - f£ e(x)dx = -iJS - 2(g + l)E(2y/g/(g + 
1)) is the ground state energy per spin. Besides, one 
can show that j g = 7r(l + As known, the total 

magnetization per spin (cr^) can be served as the order 
parameter for Ising model in a transverse magnetic field 
% (3. This supports the statement @, H S, Q that the 
geometric phase can be treated as the order parameter 
for QPT. 

In Figs. [21 [3] the real and imaginary part of the overall 
geometric phase and its derivative as functions of external 
magnetic field ft and decay parameter S are depicted. As 
can be observed j g is a continuous function of ft, if S = 0, 
and it behaves as a step-like function, if S > 0. In the 
limit cases \g\ <C 1 and \g\ > 1 we have Re7 g — > tt and 
Re7 g — > 0, respectively. 

In according to the Ehrenfest classification, the QPT 
occurred at the exceptional point, which actually is the 
circle ft;? + = 1, is of the first order QPT. In absence 
of dissipation (S = 0), we have the second order QPT. 
Indeed, as can be observed in Figs. [2j [3l the first deriva- 
tive of the ground energy (or, equivalent ly, the geometric 
phase) is the continuous function of ft and its second 
derivative diverges at the critical point ft c = 1 (5 = 0). 

In summary, we establish connection between geomet- 
ric phase and QPT in generic dissipative system and 
found the relation between the geometric phase and 
ground state energy. We show that the critical point 
where QPT occurs can be identified as the degeneracy 
point in the parameter space. Studying the critical 
behavior of the dissipative one-dimensional Ising chain 
in a transverse magnetic field, we find that the related 
QPT is of the first order QPT. In absence of dissipation 
it becomes the second order QPT. Our results support 
the claim that the relation between QPTs and geometric 
phase is a very general result, and the geometric phase 
may be considered as a good candidate to an universal 
order parameter for quantum phase transitions [4|, [5[ . 
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